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The quasi-bound modes localized on stable periodic ray orbits of dielectric micro-cavities are 
constructed in the short-wavelength limit using the parabolic equation method. These modes are 
shown to coexist with irregularly spaced "chaotic" modes for the generic case. The wavevector 
quantization rule for the quasi-bound modes is derived and given a simple physical interpretation in 
. terms of Fresnel reflection; quasi-bound modes are explictly constructed and compared to numerical 

results. The effect of discrete symmetries of the resonator is analyzed and shown to give rise to 
) quasi-degenerate multiplets; the average splitting of these multiplets is calculated by methods from 

. quantum chaos theory. 
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There is currently a great deal of interest in dielectric micro-cavities which can serve as high-Q resonators by 
q . confining light on the basis of multiple reflections from a boundary of dielectric mismatch jl], |[ || M. Such 
resonators have been used to study fundamental optical physics such as cavity quantum electrodynamics || Q and 
have been proposed and demonstrated as the basis of both active and passive optical components || ^| . Of particular 
■ interest in this work are dielectric micro-cavity lasers which have already been demonstrated for a wide variety 
i Oh) of shapes: spher esjlfjj , cylinders fjjl, squares |i"ll , hexagons and deformed cylinders and spheres (asymmetric 
resonant cavities ||~^ARCs) [|[ |T3["|l4|, |l5| , |l6| , . The work on ARC micro-cavity lasers has shown the possibility of 
producing high power directional emission from such lasers, which lase in different spatial mode patterns depending 
on the index of refraction and precise shape of the boundary. For example modes based on stable periodic ray orbits 
with bow-tie Q and triangular geometry |18 have been observed, as well as whispering gallery-like modes [[l(i| [l^] 



which are not obviously related to any periodic ray orbit. In addition a periodic orbit mode can be selected for lasing 
even if it is unstable [|l7| |l^, [2(| ; due to the analogy to quantum wavefunctions based on unstable classical periodic 
orbits, such modes have been termed "scarred" pi] . There is at present no quantitative understanding of the mode 
^v^j ■ selection mechanism in ARCs and theory has tended to work backwards from experimental observations; however 
the passive cavity solutions have been found to explain quite well the observed lasing emission patterns. The formal 
*"^» . analogy between ARCs and the problem of classical and quantum billiards has given much insight into their emission 
properties pi [|. For example it was predicted and recently observed that polymer ARCs of elliptical shape (index 
n = 1.49) have dramatically different emission patterns from quadrupolar shaped ARCs with the same major to minor 
axis ratio |^2| . The difference can be fully understood by the different structure of the phase space for ray motion 
) in the two cases, the ellipse giving rise to integrable ray motion and the quadrupole to partially chaotic (mixed) 
Oh' dynamics. Resonators are of course open systems for which radiation can leak out to infinity. In the discussion of 
the ray-wave correspondence immediately following we neglect this leakage assuming only perfect specular reflection 
of light rays at the dielectric boundary. However the theory of ARCs || || includes these effects and they will be 
' treated when relevant below. 

As is well-known, the geometric optics of a uniform dielectric region with perfectly reflecting boundaries is formally 
analogous to the problem of a point mass moving in a billiard and hence we may use the terminology of Hamiltonian 
dynamics to describe ray motion in the resonator. As discussed below, we shall specialize to cylindrical geometries in 
which the relevant motion will be in the the plane transverse to the axis, hence we focus on two-dimensional billiards. 
The quadrupole billiard is an example of a generic deformation of the circle or cylinder in that it is smooth and 
analytic and does not preserve any constant of motion. (Motion in the circle of course conserves angular momentum 
and motion in a closed elliptical cavity conserves a generalized angular momentum, the product of the instantaneous 
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angular momenta with respect to each focus (23|). A generic deformation such as the quadrupole will lead to a 
phase space for ray motion which has three types of possible motion (depending on the choice of initial conditions): 
oscillatory motion in the vicinity of a stable periodic ray orbit, chaotic motion in regions associated with unstable 
periodic ray orbits, and marginally stable motion associated with families of quasi-periodic orbits (motion on a so- 
called KAM torus). To elucidate the structure of phase space it is conventional to plot a number of representative 
trajectories in a two-dimensional cut through phase space, called the surface of section (Fig. [I]). For our system the 
surface of section corresponds to the boundary of the billiard and the coordinates of the ray are the polar angle <fi 
and the angle of incidence sin% at each bounce. The three types of motion described above are illustrated for the 
quadrupole billiard in Fig. [H both in phase space and in real space. 




FIG. 1: Surface of section illustrating the different regions of phase space for a closed quadrupole billiard with boundary given 
by r((f>) = R(l + e cos 2(f)) for e = 0.072. Real-space ray trajectories corresponding to each region are indicated at left: a) A 
quasi-periodic, marginally stable orbit, b) A stable four-bounce "diamond" periodic orbit (surrounded by stability "islands" in 
the SOS) c) A chaotic ray trajectory. Orbits of type (b) have associated with them regular gaussian solutions as we will show 
below. 

It has been shown p4j that the solutions of the wave equation for a generic shape such as the quadrupole can be 
classified by their association with these three different kinds of motion. The ray-mode (or wave-particle) correspon- 
dence becomes stronger as we approach the short-wavelength (semi-classical) limit, which in this work is defined by 
kl 3> 1 where k is the wavevector and / is a typical linear dimension of the resonator, e.g. the average radius. The 
modes associated with quasi-periodic families can be treated semiclassically by eikonal methods of the type intro- 
duced, e.g. by Keller |2f|, and referred to in its most general form as EBK (Einstein-Brillouin-Keller) quantization. 
The individual modes associated with unstable periodic orbits and chaotic motion cannot be treated by any current 
analytic methods (although the density of states for a chaotic system can be found by a sophisticated analytic method 
based on Gutzwiller's Trace Formula j26j). Finally, the modes associated with stable periodic orbits can be treated 
by generalizations of gaussian optics and will be the focus of the current work. 

If a mode is found by numerical solution, its interpretation in terms of the ray phase space can be determined with 
reasonable accuracy by means of the Husimi projection onto the phase space (see Fig. ^), although it is well-known 
that for kl not much greater than unity the exact solutions tend to smear out in the phase space over regions of order 
1/kl and do not correspond very closely to specific classical structures. For a closed generic ARC the full spectrum will 
look highly irregular (see Fig. 0(a)), but contained in the full spectrum will be regular sequences associated with tori 
and stable periodic orbits (FigTg(b)). The stable periodic orbit modes will give the simplest such sequences consisting 
of two different constant spacings, one associated with the longitudinal quantization of the orbit (free spectral range) 
and the other associated with transverse excitations. In the example of Fig. |^ the imbedded regular spectrum is due 
to the stable "bow-tie" orbit. The regular portion of the spectrum is extracted by weighting each level by the overlap 
of its Husimi function with the islands corresponding to the stable periodic orbit in the surface of section. Clearly, 
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hidden within this complex spectrum are simple regular mode sequences of the type familiar from Gaussian optics. 
In the current work we show how to calculate the resonant energies and spatial intensity patterns of such modes 
associated with arbitrary stable periodic ray orbits for both the ideal closed resonator and a dielectric resonator of 
the same shape with arbitrary dielectric mismatch n. We shall refer to these as periodic orbit modes or PO modes. 




9 

FIG. 2: Black background gives the surface of section for the quadrupole at e = 0.17 for which the four small islands correspond 
to a stable bow-tie shaped orbit (inset). A numerical solution of the Helmholtz equation for this resonator can be projected 
onto this surface of section via the Husimi transform [27] and is found to have high intensity (in false color scale) precisely on 
these islands, indicating that this is a mode associated with the bow-tie orbit. 

In the case of the open resonator, the modes have a width which can be expressed as a negative imaginary part 
of k, and some of the PO modes may be so broad (short-lived) that they would not appear as sharp spectral lines. 
Within the gaussian-optical theory we present below that this width is entirely determined by Fresnel reflection at 
the interface and would be zero for a periodic orbit which has all bounces above the total internal reflection condition, 
but would be quite large for a periodic orbit, such as the two-bounce Fabry-Perot orbit, which has normal incidence 
on the boundary. An exact solution must find a non-zero width for all PO modes, due to evanescent leakage across a 
curved interface, even if all the bounces satisfy the total internal reflection condition. 

Another limitation of the gaussian theory of stable PO modes is that it predicts exactly degenerate modes when 
the associated orbit has discrete symmetries, even in cases for which a group-theoretic analysis shows that there can 
be no exact symmetries (this is the case, for example in the quadrupole). Instead the exact solutions will have some 
integer quasi- degeneracy in which the spectrum consists of nearly degenerate multiplets, whose multiplicity depends 
in detail on the particular PO mode. This point is illustrated by the inset to Fig. ||(b). We will show below how to 
calculate the multiplicity of these quasi-degeneracies for a given PO and introduce a theoretical approach to estimate 
the size of the associated splittings. 
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FIG. 3: (a) Vertical lines indicate wavevectors of bound states of the closed quadrupole resonator for e = 0.17; no regular 
spacings are visible, (b) Spectrum weighted by overlap of the Husimi function of the solution with the bow-tie island as in 
Fig. |^. Note the emergence of regularly spaced levels with two main spacings Aki ong and Afc trans . These spac ings, indicated by 
the arrows, are calculated from the length of the bow-tie orbit and the associated Floquet phase (see section II D below). The 
color coding corresponds to the four possible symmetry types of the solutions ( see se ction |v| below). In the inset is a magnified 
view showing the splitting of quasi-dege nerat e doublets as discussed in section [VC. Note the pairing of the (+, +) and (+, — ) 
symmetry types as predicted in section [VB. The different symmetry pairs alternate every free spectral range (Aki ong ). 



II. GAUSSIAN OPTICAL APPROACH TO THE CLOSED CAVITY 



The quantization of electromagnetic modes within dielectric bodies enclosed by metallic boundaries, including the 
case of arbitrary index variation inside the resonator and three-dimensions has been treated before in the literature 
ps| . However the generality of the treatment makes it difficult to extract simple results of use to researchers working 
with uniform dielectric optical micro-cavities; more importantly all the work of which we are aware focuses on the case 
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of Dirichlet boundary conditions corresponding to perfect reflection at the boundary. Perfect reflection of course leads 
to true bound states. In the next section we show how to generalize these results to the correct boundary conditions 
at a dielectric interface and hence for the case of quasi-bound as opposed to bound states. However, first, in this 
section we develop the formalism for the closed case which we will generalize to the case of interest. In all of this 
work we will specialize to the case of two dimensions, corresponding to an infinite dielectric cylinder with an arbitrary 
cross-section in the transverse plane (see Fig. |j) with the condition k z — 0. To conform with conventions introduced 
below we will refer to the two-dimensional coordinate system as (X, Z) . In this case the TM and TE polarizations 
separate and we have a scalar wave equation with simple continuity conditions at the boundary for the electric field 
(for TM) or magnetic field (for TE). For the closed case (for which the field is zero on the boundary) we can set the 
index n = 1 and work simply with the Helmholtz equation. 




FIG. 4: Illustration of the reduction of the Maxwell equation for an infinite dielectric cylinder to the 2D Helmholtz equation 
for the TM case (E field parallel to axis) and k z — 0. 

Consider the solutions E{X, Z) of the Helmholtz equation in two dimensions: 

(V 2 + fc 2 )£ = (1) 

This solution is assumed to be defined in a bounded two dimensional domain D, with a boundary dD on which E = 0, 
leading to discrete real eigenvalues k. We are interested in a subset of these solutions for asymptotically large values of 
k for which the eigenfunctions are localized around stable periodic orbits(POs) of the specularly reflecting boundary 
dD. 



A. The Parabolic equation approximation 



The "N-bounce PC's corresponding to a boundary dD are the set of ray orbits which close upon themselves upon 
reflecting specularly N times. The shape of the boundary defines a non-linear map from the incident angle and polar 
angle at the m th bounce (</> m , sin^m) to that at the m + I bounce. Typical trajectories of this map are shown in the 
surface of section plot of Fig. 1. The period-N orbits are the fixed points of the N th iteration of this map. For a given 
period-N orbit (such as the period four "diamond" orbit shown in Fig. ||), let the length of mth segment ("arm") 
be l m , the accumulated distance from origin be Lm — J2m=i 'mi an d L = Ln be the length of the entire PO. We 
are looking for modal solutions which are localized around the PO and decay in the transverse direction, hence we 
express Eq. (Q) in Cartesian coordinates (x m ,z m ) attached to the PO, where z m -axis is aligned with mth arm and 
x m is the transverse coordinate. We also use z to denote the cumulative length along the PO, which varies in the 
interval (— oo,+oo). 

We write the general solution as: 



N 



E(X, Z)=Y, E m (x m (X, Z), z m (X, Z)) (2) 

m=l 

where the "local set" (x m ,z m ) and the fixed set (X,Z) are related by shifts and rotations (see Fig. ||). Next, in 
accordance with the parabolic equation approximation |2g| ], we assume that the main variation of the phase in z- 
direction is linear ("slowly varying envelope approximation") and factor it out: 

Emi^miZm) — ^m(^m) ^m)^ (3) 

This definition suggests defining the origins of the local coordinates such that z m = z along the PO, and we do so 
(see Fig. |). 

Next, we insert the solution Eq. (|^) into Eq. (Q) and using the invariance of the Laplacian, we obtain: 

(V 2 m + k 2 ) E m = (4) 
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FIG. 5: Coordinate system and variables used in the text displayed for the case of a quadrupolar boundary dD and the diamond 
four-bounce PO. A fixed coordinate system (X,Z) is attached to the origin. The "mobile" coordinate systems (x m ,z m ) are 
fixed on segments of the periodic orbit so that their respective z-axes are parallel to the segment, while their origins are set back 
a distance L m (or nL m for transmitted beam axes), so as to account for zeroth order phase accumulation between successive 
bounce-points. £1,^2 are the common local coordinates at each bounce (index m suppressed). Scaled coordinates are denoted 
by tildes, e.g. x m = \fkx m . The coordinate transformations at each bounce m are given by Zi = L m + £1 sinxi + £2 cosxii 
z r = L m + £isinxi - ^cosxi, Zt = nL m + £isinxt + ^cosxt and Xi = £icosXi - £2shiXi, x T = £icosxi + ^sinxi, 
xt — £1 cos xt — £2 sin xt, where i, r, t refer to the incident, transmitted and reflected solutions. 



and the boundary condition translates into 



En 



E, 



m+l\dD 







(5) 



Here V m is the Laplacian expressed in local coordinate system. This reduction is possible as long as the solutions 
are well-localized, and the bounce points of the PO are well-separated (with respect to even if the PO were to 

self-intersect. These assumptions will be justified by the ensuing construction. 

Dropping for the moment the arm index, we will focus on Eq. (|). Inserting Eq. (|), we arrive at 



2iku, 







(6) 



The basic assumption of the method is that after removing phase factor exp[ifcz] which varies on the scale of the 
wavelength A, the z-dependence of u is slow, i.e. u z ~ u/l, u zz ~ u/l 2 , where I is a typical linear dimension associated 
with the boundary, e.g. a chord length of the orbit or the curvature at a bounce point. In the scmiclassical limit / ^> A. 
The transverse (x) variation of u on the other hand is assumed to occur on a scale vTX, intermediate between the 
wavelength and the cavity scale, I] hence the transverse variation of u(x, z) is much more rapid than its longitudinal 
variation. This motivates us to introduce the scaling x = ykx to treat this boundary layer, which leads to 



mx + \u zz + 2iu z = 0. 
k 



(7) 



We can now neglect the u zz term due to the condition kl 1, and obtain a partial differential equation of parabolic 
type: 



C u{x, 







(8) 



where C = 9| + 2id z . Next, we make the ansatz 



u{x, z) = c A(z) exp 



-n(z)x 2 
2 v ; 



(9) 
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Inserting Eq. (^) into Eq. (||) and requiring it to be satisfied for all x, we obtain the relations 

fl 2 + Q' = (10) 
An + 2A' = (11) 

Here and in the rest of the text we will use primes as a shorthand notation for a z-derivative. Next, making the 
substitution 

Q{z) 

we obtain 

Q" = (13) 
n' a 1 

| + 2-^„ ,14, 

Note that Eq. ( |l3| ) is the Euler equation for ray propagation in a homogeneous medium with general solution Q(z) = 
az + (3; we will be able to interpret Q(z) below as describing a ray nearby the periodic orbit with relative angle and 
intercept determined by a, /?. 

B. Boundary Conditions 

Having found the general solution of Eq. (||) along one segment of the periodic orbit we must impose the boundary 
condition Eq. (^J) in order to connect solutions in successive segments. Writing out this condition: 



7= , f exp ikz m + -Cl(z m )x m + == exp ikz m+l + -Cl(z m+ i)x m+1 \ dD 

\/Qm\Z m ) \ 1 J \/Qm+l{z m +lj X 



This equation must be satisfied on an arc of the boundary of length ~ yAZ around the reflection point. Since this 
length is much smaller than I we can express this arc on the boundary dD as an arc of a circle of radius p (the 
curvature at the reflection point). We express the boundary condition in a (scaled) common local coordinate system 
for the incident and reflected fields (£i, £2) pointing along the tangent and the normal at the bounce point (see Fig. |jj). 
Because the boundary condition must be satisfied on the entire arc it follows the phases of each term must be equal, 

k (l m + -^=|i sinx 




k U m + — £1 sin Xm + TT~ rus X'" 1 ~ vV" " ' j ros \-« + T.~ — sin \<- 
and there is a amplitude condition as well, 



= 



\J Qm(L m ) y Qm+l{L m ) 

Here, we assume that = O(l), and we will carry out the solution of these equations to 0(1/ yk). Note that, it 
is sufficient to take Q{z)\qd ~ Q(lm)j a t this level. In each segment we have three constants which determine our 
solution: c m , a m , (3 m (where Q m (z) — a m z + (3); however due to its form our solution is uniquely determined by the 
two ratios f3 m /a m , c m j ^Ja m . Therefore we have the freedom to fix one matching relation for the Q m by convention, 
which then determines the other two uniquely. To conform with standard definitions we fix Q m +i = Q m , then the 
amplitude and phase equality conditions give the relations: 

Q'm+l J V ~ o„, cos v™ 1 / V Q'm J ™ V Qm J ^ 



and 



C m+1 = e-™C m (16) 
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Here we always choose the principal branch of \fQ. These conditions allow us to propagate any solution through a 
reflection point via the reflection matrix lZ m and Eq. (|l6| ) . Note that with these conventions the reflection matrix is 
precisely the standard "ABCD" matrix of ray optics for reflection at arbitrary incidence angle (in the tangent plane) 
from a curved mirror |p9| . It follows that Q m and Q m +i can be interpreted as the transverse coordinates of the 
incident and reflected rays with respect to the PO, and that then Eq. ( |l5| ) describes specular ray reflection at the 
boundary, if the non- linear dynamics is linearized around the reflection point. 



C. Ray dynamics in phase space 

To formalize the relation to ray propagation in the paraxial limit we define the ray position coordinate to be Q and 

the conjugate momentum P = Q', with z playing the role of the time. Let's introduce the column-vector (p). The 
"fundamental matrix" II is the matrix obtained by two such linearly independent vectors 

where the linear independence, as expressed by the Wronskian condition W(Q', Q) 7^ reduces to detll 7^ 0. Then, 
for example, the Euler equation Eq. ([l3]) in each arm can be expressed as 

^ = jnn (is) 

dz 

where H = (00) an( A J = ( 0-1 ) ■ It's straightforward to show that detll is a constant of motion. To take into account 
the discreteness of the dynamics in a natural manner, we will introduce "coordinates" for the ray, represented by a 
z-independent column vector h. Any ray of the mth arm in the solution space of Eq. (Oh can then be expressed by 
a z-independent column vector h m 



Q{z) 

P(z) 



= n(z)h m (19) 



We will choose H(z) = (01), so that II 1 (z) = n(— z). Then if Q m (z) = a m z + /3 m , we have h m = ( a™ )■ In this 
notation, propagation within each arm (i.e. z + I < L m where L m -i < z < L m ) is induced by the 'propagator' II 

^ + l)\- m (%^)\ (20) 



Pm{z + l) J V > \ Pm{z) 

which is the "ABCD" matrix for free propagation p9[ . Thus, for example z —> z' propagation for L TO _i < z < L rn < 
z' < L m+ i is 

Qm+l («') ^ = n(z , _ Lm)UmU{Lm _ z) ( Qm{z) \ (21) 



Of special importance is the monodromy matrix, A4(z), which propagates rays a full round-trip, i.e. by the length L 
of the corresponding PO and is given by 

M{z) = n(z - L m -i)n m - l n.{i m -i) ■ ■ ■ ii(i m+1 )K m TL(L m - z) (22) 

for £ m -i < z < L m . Note that l m +N — lm and 1Z m +N = Although the specific form of M.(z) depends on the 
choice of origin, z, it is easily shown that all other choices would give a similar matrix and hence the eigenvalues of 
the monodromy matrix are independent of this choice. We will suppress the argument z below. 



D. Single-valuedness and quantization 

Having determined how to propagate an initial solution Q, P an arbitrary distance around the periodic orbit we 
can generate a solution of the parabolic equation which satisfies the boundary conditions by an arbitrary initial choice 
Q(0),P(0). However an arbitrary solution of this type will not reproduce itself after propagation by L (one loop 
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around the PO). Recalling that our solution is translated back into the two-dimensional space (X, Z) by Eq. (0), and 
that the function E(X, Z) must be single- valued, we must require periodicity for our solutions 

E(x,z + L) =E(x,z). (23) 

We will suppress again the reference to arm index and use the notation E(x,z) — E m (x m , z m ) and u(x, z) = 
u m (x m , z m ) whenever L m _i < z < L rn . Since the phase factor in E advances by exp[i£;L] with each loop around the 
periodic orbit and single-valuedness implies that 

u(x,z + L)e ikL = u(x,z). (24) 

This periodicity condition will only be solvable for discrete values of k and will lead to our quantization rule for the 
PO modes. 

From the form of u in Eq. ([)]) we see that the phase £l(z) = Q'/Q will be unchanged if we choose (Q,P) to be an 
eigenvector (<?i,pi)), (92,^2) of the monodromy matrix as in this case the Q'(z + L) = Xx^Q 1 [z) , Q{z + L) = \x^Q{z) 
and the ratio £l(z) = £l(z + L). The monodromy matrix is unimodular and symplectic and its eigenvalues come in 
inverse pairs, which are either purely imaginary (stable case) or purely real (unstable and marginally stable cases). 
If the PO is unstable and the eigenvalues are real, then the eigenvectors and hence are real. But a purely real 
Q(z) means that the gaussian factor in u(x, z) is purely imaginary and the solution does not decay in the direction 
transverse to the PO, contradicting the initial assumption of the parabolic approximation to the Hclmholtz equation. 
Hence for unstable POs our construction is inconsistent and we cannot find a solution of this form localized near the 
PO. 

On the other hand, for the stable case the eigenvalues come in complex conjugate pairs and the reality of the 
monodromy matrix then implies that the eigenvectors cannot be purely real and are related by complex conjugation. 
Therefore we can always construct a gaussian solution based on the eigenvector with ImVl > so that there is a 
negative real part of the gaussian exponent and the resulting solution decays away from the PO. It is easy to check 
that the range of the decay is ~ y/M as assumed above. 

The eigenvectors of the monodromy matrix are sometimes referred to as Floquet rays of the corresponding Euler 
equation. In a more formal discussion which can be straightforwardly generalized to arbitrary index variation and 
three dimensions psfl our condition on the choice of Q can be expressed by recalling that the two eigenvectors must 
generate a constant Wronskian, and we will choose the eigenvector [q{z),p{z)) such that 

pq* -qp* =i. (25) 

This condition implies /mO > 0. Note that the overall magnitude of the Wronskian is determined by our choice of 
normalization of the eigenvectors and is here chosen to be unity. 

Having made this uniquely determined choice (up to a scale factor) and using (p( 2 +L) ) = e llp {p(z] ) an d c m +jv = 
e _MrJV c m we obtain 

u(x, z + L) = e-^-^^+^uix, z) (26) 

From Eq. (^3|) and Eq. (^) , we finally obtain the quantization rule for the wavevectors of the bound states of the 
closed cavity: 

kL= ^-<p + 2nm + mod 2w [(iV + N^w] (27) 

where m is an integer, TV is the number of bounces in the PO, (p is the Floquet phase obtained from the eigenvalues 
of the monodromy matrix and is an integer known as the Maslov index in the non-linear dynamics literature [ |§0| . 
It arises in the following manner. As already noted, the eigenvalues of the monodromy matrix for stable POs are 
complex conjugate numbers of modulus unity whose phase is the Floquet phase. We can define <p = Aig[y/ q(L) /q(0)], 
where Arg[-] denotes the principal argument; hence the Floquet phase depends only on M. However our solution, 
Eq. (Q), involves the square root of q(z) and will be sensitive to the number of times the phase of q{z) wraps around 
the origin as z goes from zero to L. If this winding number (or Maslov index) is called then the actual phase 
advance along the PO is <p + if Nu is odd this leads to an observable it phase shift in the solution not included 

by simply diagonalizing M. to find the Floquet phase. 



Nfj_ may be directly calculated by propagating q(z): 



1 



L 

d(ln q{z)) 







(28) 
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FIG. 6: Intensity of the TM solution for a bow-tie mode plotted in a false color scale, (a) calculated numerically and (b) from 
the gaussian optical theory with parameters m = 100, <p = 2.11391, iV M = 1 and N = 4. Note the excellent agreement of the 
quantized values for kR (R is the average radius of the quadrupole). 



where [•] denotes the integer part. There is no simple rule for reading off the Maslov index from the geometry of the 
PO; however the Maslov index doesn't affect the free spectral range or the transverse mode spacing. 

Except for the Maslov index, the quantization rule Eq. (|2^) is familiar from Fabry-Perot resonators. The longitudinal 
mode index m gives rise to a free spectral range Afc; on9 = 2ir/L for gaussian modes of a stable PO of length L. The 
Floquet phase (p/2L is the zero-point energy associated with the transverse quantization of the mode and we will 
shortly derive excited transverse modes with spacing ip/L. In Fig. ^, we plot for comparison the analytic solutions 
for the bow-tie resonance just derived in comparison to a numerical solution of the same problem; both the intensity 
patterns and the quantized value of k agree extremely well. 



E. Transverse excited modes 

As is well known, the solution Eq. ^ is only one family of the possible solutions of the parabolic equation Eq. (H) 
satisfying the boundary conditions Eq. (^|) and the periodicity condition Eq. (p3|), the one corresponding to the ground 
state for transverse oscillations. Further solutions can be generated with algebraic techniques which were originally 
developed in the context of quantum oscillators Q . If we refer to Eq. (||) we see that we are looking for solutions of 
C(u) — 0, i.e. eigenfunctions of the differential operator with eigenvalue zero. It is natural following the analogy to 
quantum oscillators to seek additional solutions by defining lowering and raising operators. 

A(z) = -iq(z)d s +p(z)x (29) 

A*(z) = iq*{z)d i +p*(z)x. (30) 

We can easily show that the operators (A, A',£) form an algebra. Namely, that [At,£] = [A, C] = and furthermore 
that [A, A^] = i(p*q—q*p) = 1. Defining the "ground-state" solution we have found as u^°\ the commutation condition 
implies that (A^u^)(x, z) is also a solution of C(u) = 0. Further it can be checked that while (Au^) = 0, (A^)u^ 
is a non-trivial solution. We can find the wavevector quantization condition for this solution by calculating the 
additional phase acquired by (A^u^)(x, z) upon performing one loop around the orbit (A^u^°')(x, z + L). Noting that 
q*(z +L) = e- llfi q*(z),p*(z + L) = e~ ilfi p*(z), we find that A^(z + L) = e'^A^z). Thus this solution will acquire an 
additional phase —tp with respect to This means that the I th family of solutions 



U^(x,z) = (At)V°) 



(31) 
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will satisfy the wavevector quantization rule: 



kL=(l + -)tp + 2irm + mod 2 ,r[(iV + N^ic) 



(32) 



This result has the well-known interpretation of adding I transverse quanta kx = f/L to the energy of the ground 
state gaussian solution. Thus (in two dimensions) the general gaussian modes have two modes indices (l,m) corre- 
sponding to the number of transverse and longitudinal quanta respectively and two different uniform spacings in the 
spectra (see Fig. ||(b)). 

In order to obtain explicit forms for the excited state solutions found here one may use the theory of orthogonal 
polynomials to find 




(x,z) 



(33) 



where Hi are the Hcrmitc polynomials. 

As noted in the introduction, the solutions we have found by the parabolic equation method do not reflect correctly 
the discrete symmetries of the cavity which may be present. The theory of the symmetrized solutions and the breaking 
of degeneracy is essentially the same for both the closed and open cavity and will be presented in section (IV) after 
treating the open case in the next section. 



III. OPENING THE CAVITY - THE DIELECTRIC RESONATOR 



We now consider the wave equation for a uniform dielectric rod of arbitrary cross-section surrounded by vacuum. 
As noted in the introduction, Maxwell's equations separate into TM and TE polarized waves satisfying 

(V 2 + n(r) 2 fc 2 ) * = (34) 

Maxwell boundary conditions translate into continuity of and its normal derivative across the boundary 3D, defined 
by the discontinuity in n. Here VP = E(^ = B) for TM(TE) modes. We will only consider the case of a uniform 
dielectric in vacuum for which the index of refraction n(r 6 D) = n and n(r & D) = 1. Thus we have the Helmholtz 
equation with wave vector nk inside the dielectric and k outside. The solutions to the wave equation for this case 
cannot exist only within the dielectric, as the continuity conditions at the dielectric interface do not allow such 
solutions. On physical grounds we expect solutions at every value of the external wavevector k, corresponding to 
elastic scattering from the dielectric, but that there will be narrow intervals 5k for which these solutions will have 
relatively high intensity within the dielectric, corresponding to the scattering resonances. A standard technique for 
describing these resonances as they enter into laser theory is to impose the boundary condition that there exist only 
outgoing waves external to the cavity. This boundary condition combined with the continuity conditions cannot be 
satisfied for real wavevectors k and instead leads to discrete solutions at complex values of k, with the imaginary part 
of k giving the width of the resonance. These discrete solutions are called quasi-bound modes or quasi-normal modes 
p2[ . We now show how such quasi-bound modes can be incorporated into the gaussian optical resonance theory just 
described. 

As before we look for solutions which, within the cavity, are localized around periodic orbits based on specular 
reflection within the cavity. In order to satisfy the boundary conditions the solutions outside the cavity will be 
localized around rays extending out to infinity in the directions determined by Snell's law at the bounce point of the 
PO. The quasi-bound state condition is imposed by insisting that the solutions along those rays to infinity are only 
outgoing. This can only be achieved by making k complex (assuming real index n). It is worth noting that there are 
well-known corrections to specular reflection and refraction at a dielectric interface, for example the Goos-Hanchen 
shift [ p3| . It is interesting to attempt to incorporate such effects into our approach at higher order, however we do 
not do so here and confine ourselves to obtaining a consistent solution to lowest order in kl. 

As before we will define the solution as the sum of solutions E m (x m , z m ) attached to each segment of the PO and 
define local coordinates (x m , z m ) attached to the PO. Now in addition we need an outside solution at each reflection 
point E mt with its own coordinate system (x m t , z m t) rotated by an angle given by Snell's law applied to the direction 
of the incident ray (see Fig. ||). The ansatz of Eq. (^) thus applies, where now the sum will run over 2N components, 
which will include the transmitted fields. Introducing the slowly varying envelope approximation Eq. (||) and the 
scalings x mt = V~kx, x m = Vnkx, we get the parabolic equation Eq. (^) for each component, at lowest order in k. 
The boundary conditions close to the m bounce point will take the form 



Ei + E r \ dD E t \g D 



(35) 
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and 



d n Ei + d n E r \g D - — d n Et\g]j+ 



(36) 



Here d n is the normal derivative at the boundary. The alternative indices i, r, t stand for m, m+1 and mt, respectively. 
Since the parabolic equation Eq. (0) is satisfied in appropriately scaled coordinates within each segment, we write all 



solutions in the general form Em = Am exp(i<$>M ) where $ 



M 



riMkz - 



Qjx 2 M and A 



si = c m /VQM- Here M 
and cm, so that the 



' M 



stands for m or mt and n m — n, n mt — 1. As for the closed case, we need to determine Qm 
boundary conditions are satisfied, and then impose single-valucdness to quantize k. 

Similarly to the closed case, the first continuity condition Eq. (|35| ) must be satisfied on an arc of size ~ y/Xl on the 
boundary around each bounce point and that implies that the phases of the incident, transmitted and reflected waves 
must be equal. This equality will be implemented in the coordinate system (£1,^2) along the tangent and normal to 
the boundary at the reflection point, as before. We can again expand the boundary as the arc of a circle of radius /?, 
the curvature at the bounce point. Since the equations are of the same form for each reflection point it is convenient 
at this point to suppress the index m and use the indices i,r,t to denote the quantities associated with the incident, 
reflected and transmitted wave at the m th bounce point. 

The analysis of the phase equality on the boundary for the incident and reflected waves is exactly the same as in 
the closed case and again leads to Eq. ( |l5| ) describing specular reflection. Equating the phases of the incident and 
transmitted wave leads to: 



® t = k\ nL r . 



1 ? . 1 CI \ , 1 Q'i 



1 i ■ 1 Ci 2 

=£1 sin xt r 7T cos Xt 

Ink nk 2p 



fiCOSx+ -y^sinx 
nk Zp 



1 Q't h a. 1 ^ 2 • 

2^UiCOS Xt + --sm Xt 



where we have defined x — Xmi P = Pm and all quantities are evaluated at z = L m . Recalling that nsmxi = sin^t 
we get up to 0(l/vfe) 



Qi 
Q't 



i/n 



2(1-M) 



(37) 



where fi — cos\i/ cosxt and the relation Qi — pQ t is a convention similar to Q r — Qi. Again, the matrix in Eq. ( |3^ ) 
is just the ABCD matrix for transmission of rays through a curved dielectric interface at arbitrary angle of incidence 
in the tangential plane |^9| . 

Using the phase equality on the boundary the continuity of the field gives the general transport equation: 



(-1 



\Qi \Qr V 

which becomes (using the conventions Q r — Qi , Qt — Qi / p) 

Ci + C r = y/flCt 



Qt 



(38) 



(39) 



In order to find the quantization condition we need a direct relation between Cj and c r as we had in the closed case. 
This is provided by the normal derivative boundary condition Eq. (|3^). Keeping only the leading terms this condition 
becomes: 



C 4 9„$i + c r <9„$ r = yfpCtdn^t 



(40) 



At the level of approximation needed one finds the simple results <9„<I>.; = nfccosx, d n $ r = —nkco&x,d n ^ t = kcosxt, 
leading to 



and 



ny/jl(Ci - c r ) — c t 
np — 1 



np + 1 



-Ci- 



(41) 



(42) 
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Note the key result that 



M 2 



\ncosxi - cosxt| s 
\ncos\i + cosxtp 



M 5 



which is precisely the Fresnel reflection law at a (flat) dielectric interface. 

Now we impose the single- valuedness or periodicity condition to obtain the quantization rule for k. 

E(x,z + L) =E(x,z) 



(43) 



(44) 



Note however that we have a qualitatively different situation than for the closed cavity; some amplitude can be lost 
at each reflection and it will in general be impossible to make a loop around the PO and return to the same field 
amplitude unless k is complex. We have the condition 



u(x, z + L)e inkL 



u(x, z) 



(45) 



which can be satisfied by choosing Q(z), P(z) to be the appropriate eigenvector of the monodromy matrix (note that 
M. is unchanged from the closed case as it only pertains to the propagation of the phase) and with this choice the 
quantization condition becomes 



1 



JY 



ikL — —(f + 27rm + mod27r[(^V + N^ir] — i log 

6=1 



n^b - 1 



n(j,b 



1 



(46) 



Recall that the Fresnel reflection law has the property that it gives a pure phase for rays incident above total 
internal reflection. Thus the new term in the quantization law due to Fresnel reflection can be either purely real (all 
bounces of PO totally internally reflected), purely positive imaginary (all bounces below T1R) or complex (some TIR 

then 



bounces, some refracted bounces). If we define: iff — Re[— i^^log 
the quantization rule gives 



nfi b -l 



] and 7/ = Im[-i £f log 



Re[nkL] 
ImlnkL] 



2irm + mod 27r [( N + N^tt] + tp/2 + ip f 
-If- 



(47) 
(48) 



As noted above, this result is only in the leading order approximation, and it ignores both the effects of evanescent 
leakage at a curved interface and the momentum width of the gaussian "beam" which leads to violations of ray optics. 
These effects will give a non-zero imaginary part (width) to all resonances, even those with all bounce points above 
TIR. 

In table | we present a comparison between the numerically obtained quantized wavevectors for a bow-tie resonance 
and the values for the real and imaginary part of k predicted by Eq. ([h|) for three different indices of refraction. Note 
that the best agreement is for the case far from total internal reflection, and the worst agreement is the case near 
TIR. 



index 
of refraction 



numerical 
calculation 



gaussian 
quantization rule 



surface of 
section 



n = 2.0 



n = 2.9 



5.1 



nkR = 100.53788 + 0.49758i 



nkR = 100.53858 + 0.49635i 



nkR = 100.59376 + 0.16185i 



nkR = 100.53858 + 0.14178i 



nkR = 100.84716 + 0.00245i 



nkR = 100.85111 + O.OOOOOi 



— 



pin 



TABLE I: Table comparing the gaussian optical prediction for the complex k values of bow-tie resonances of the quadrupole 
(e = 0.17) to numerically obtained values for three different indices of refraction corresponding to incidence below, at and above 
the critical angle for total internal reflection. The schematic at right depicts the bow-tie island in comparison to the critical 
line for total internal reflection for the three cases (horizontal lines). 
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IV. SYMMETRY ANALYSIS AND QUASI-DEGENERACY 



As noted above, the gaussian theory we have just presented predicts exact degeneracies if there exist several 
symmetry-related stable orbits (note that in this context the same path traversed in the opposite sense is considered a 
distinct symmetry-related orbit). However it is well-known that wave equations with discrete symmetries cannot have 
a degeneracy which is larger than the largest dimension of the irreducible representations of the symmetry group of 
the equation. Here we are concerned with the point group (rotations and reflections) of the dielectric resonator in two 
dimensions. Let G be a group which leaves the cavity invariant. Then, if ^(x) is a solution to the Helmholtz equation, 
so is 'l'(gx), where g £ G and x = (X,Z). The symmetrized solutions, i.e. solutions which transform according to 
the irreducible representations of G under the action of G, can be obtained by the projection operators of the group: 

V m E(x) = f%J2 Xrn(9)E(gx) (49) 

Here d m is the dimensionality and Xmig) is the character of the m th irreducible representation and the solution so 
obtained, denoted by E m (jx) is the resulting symmetry-projected solution. For a given irreducible representation there 
are as many symmetrized solutions as the dimension of that irreducible representation of G. We will focus here on 
the case of the closed resonator, but the general principles apply to the open case as well. 



A. Symmetrized modes for the quadrupole 

Let's consider our canonical example, the quadrupole (see Fig. |J). The symmetry group of the quadrupole is G = 
C2® C2 — {1, <?x, cz, oxoz}-, the group of reflections about the X and Z axes. This group has four one-dimensional 
representations only, and thus cannot have any exactly degenerate solutions (barring accidental degeneracy). The 
existence of four irreducible representations means that given the one solution E(X, Z) we have constructed to the 
Helmholtz equation, we can generate four linearly independent solutions by projection according to Eq. ( |49| ) above. 
We will label the representations by m — (r, s), where r,s — ± denotes the action of inversion of X and Z respectively. 
The symmetrized solutions are then 



where 



E( ++ ) = \{ei + e 2 + e 3 + e 4 ), £(+_) = \{e x + e 2 - e 3 - e 4 ) 
E(-+) = \{e\ - e 2 + e 3 - e 4 ), E( — ) = \{e\ - e 2 - e 3 + e 4 ) 



ei =E(X,Z), e 2 = E(-X,Z), e 3 = E{X,-Z), e 4 = E(-X,-Z). 



(50) 



(51) 



In this case the symmetrized solutions are just the solutions with definite parity with respect to the symmetry axes 
of the quadrupole. 

The key point here is that the within the parabolic equation approximation these four solutions are exactly degen- 
erate, whereas our group-theoretic analysis for the exact Helmholtz equation tells us that they cannot be so, although 
they will be nearly degenerate. Moreover our original solution E{X, Z) cannot be an exact solution, as it does not 
transform as any irreducible representation of the symmetry group. A further important point is that while we can 
always construct a number of symmetrized solutions equal to the sum of the dimensions of the irreducible representa- 
tions, there is no guarantee that such a projection will yield a non-trivial solution. In fact in the case of the quadrupole 
we will show below that for each quantized value of k only two of the projected solutions are non-trivial, leading to 
quasi-degenerate doublets in the spectrum. We will present below a simple rule which allows one to calculate the 
quasi-degeneracy given the periodic orbit and the symmetry group of the resonator. 

Before discussing the general rule, we illustrate the basic procedure for the case of a bow-tie PO. Let £\ and £2 be 
the lengths of the vertical and diagonal legs of the bow-tie, so that L = 2{£\ + £2) is the total length. Then, 

• = 1 



ei = E{gx) 



VOW 



exp 



ikz + — Q(z)x 



(52) 



• 9 = a x '■ i z — * L/2 + z, x - 



('•2 



■>-z) 



^e l ^/ 2 q{z) 



exp 



L i 
ik(z + — + -Sl{z)x 2 



(53) 
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• .9 = °~z '■ z ^ t\ — z , x —> x) 

e 3 = E{gx) 

• g = <Jx<Jz '■ {z — > i/2 + — z, x —>■ —x) 

e 4 = E(gx) = eM 



i-z) 



: exp 



ik{h -z) + - z)x 2 



■. exp 



ik(~+ti -z) + - z)x 2 



(54) 



e lC e 3 



(55) 



where the phase factor £ = 1/2 (tt — ip/2 + kL). Here we use the fact that A4l = Ml/2-Ml/2 for the bow-tie orbit, 
where Ml is the monodromy matrix for the whole length L. It follows that q(z + L/2) = e %tp / 2 q(z). Note also the 
appearance of the factors e'™, which is due to the specific choice of branch-cut for yj q(z). Putting these together, 
we obtain 



s (++) 



|(er- 



cos | 
e 3 )e 1 ^ cos | 

*i . c 



eaje 2 sin ■ 



(56) 



e 3 )e 2 sin | 



The solutions are fc-dependent and must be evaluated for the quantized values of k. Referring to the quantization 
condition Eq. (p7|) we find that the phase £ = to,7t where to is the longitudinal mode index of the state. Hence 



E, 



Ei A i oc cos 



(++),^(+-) 
Ei h \,Ei_j oc sin 







2 



m = 1, 3, 5, . 
m = 0,2,4, . 



(57) 



Thus the quasi-degeneracy of the solutions is two for the bow-tie, the solutions with identical parity under az form 
the doublets (see inset Fig. ||b)), and these two parity types alternate in the spectrum every free spectral range. Note 
that while we have illustrated the analysis for for the ground state (I = 0) one finds exactly the same result for the 
I th transverse mode, with doublets paired according to the index to, independent of I. 

This illustrates a general procedure, valid for any stable PO. First, one finds by the parabolic equation method a 
non-symmetrized approximate solution E{X, Z) localized on the PO. Second, one generates the symmetrized solutions 
from knowledge of the irreducible representations of the symmetry group. Third, one evaluates these solutions for the 
quantized values of k; the non-zero solutions give one the quasi-degeneracy and the symmetry groupings (e.g. (H — h) 
with (H — ) in the above case). The same principles apply to mirror resonators with the same symmetry group. Note 
that in the case of a high symmetry resonator (or mirror arrangement) e.g a square or a hexagon, for which there 
exist two dimensional irreducible representations, exact degeneracy is possible and can be found by these methods. 



B. Simple Rule for Quasi-Degeneracy 

Although the construction just presented allows one to find the quasi-degeneracy and symmetry pairing, it is 
convenient when possible to have a simple rule to get the quasi-degeneracy and symmetry-pairing from the geometry 
of the orbit. The quasi-degeneracy is easily determined by the following rule: 

The quasi-degeneracy D is equal to the number of distinct classical periodic orbits which are related by the spatial 
symmetry group and time-reversal symmetry. 

In this rule "distinct" orbits are defined as orbits which cannot be mapped into one another by time translation. 
Therefore a self- retracing orbit such as the Fabry-Perot, two-bounce orbit, only counts as one orbit and is non- 
degenerate (see table |n|). In contrast for a circulating orbit like the diamond no translation in time will take the orbit 
into its time-reversed partner. This rule can be obtained from semiclassical methods similar to the Gutzwiller Trace 
formula fl3~4j . The density of states can be expressed by a summation over periodic orbits and their repetitions; for 
stable periodic orbits the summation over repetitions yields a delta function at the semiclassical energies (corresponding 
to the same wavevectors as we find from our quantization rule). This approach would give an alternative derivation 
of our results which is less familiar in optics than the parabolic equation method we have chosen. However the 
semiclassical method makes it clear that there will be a mode for every distinct symmetry-related PO using the 
definition we have just given (of course in this method, as in the parabolic equation method, one would predict an 
exact degeneracy instead of the quasi-degeneracy we have discussed) . 
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orbit 



symmetry- 
related 



time- 
reversal 



quasi- 
degeneracy 



symmetry- 
reduced orbit 



symmetry- 
pairing 



bowtie 



[(++),(+-)];[(—),(-+) 



diamond 



[(++),(—)];[(+-).(-+)] 



triangle 



[(++),(+-),(-+),(" 



fish 



[(++),(-+)];[(+-),(—)] 



fabry perot 



[(++)];[(+-)] 



TABLE II: Table illustrating the application of the two symmetry rules to five simples POs. The second column is the total 
number of orbits of this shape related by spatial symmetry; the third column is the total number of orbits of this shape 
generated by time-reversal symmetry. By Rule 1 the quasi-degeneracy is the product of these two numbers. Column 4 gives 
the symmetry-reduced orbit which leads via Rule 2 to the symmetry pairing indicated in column 5. 



Let us illustrate the application of this rule. The bow-tie orbit goes into itself under all the reflection symmetries 
and so spatial symmetry generates no new orbits; however time-reversal changes the sense of traversal of each leg of 
the orbit and does give a distinct orbit. Thus the predicted quasi-degeneracy is two, which we found to be correct by 
our explicit construction above. In contrast, the triangle orbit (see table |n]) has a symmetry related distinct orbit and 
a definite sense of circulation which is reversed by time-reversal, hence it should have a quasi-degeneracy 2x2 = 4 
leading to quartets instead of doublets. A few different cases of this rule are illustrated in table [□]. 

The rule we have just given tells one the quasi-degeneracy, D, but not the symmetry-pairing. For the case of 
reflection symmetries one can state a second rule which determines these pairings. First fold the PO of interest back 
into the symmetry-reduced resonator |34 (see table ||) using reflection until it completes one period in the reduced 
resonator. The symmetry-reduced resonator has boundaries which correspond to lines of reflection symmetry in the 
original problem. Anti-symmetric solutions with respect to each of these lines of symmetry correspond to Dirichlet 
boundary conditions; symmetric solutions must have zero derivative corresponding to Neumann boundary conditions. 
The boundary conditions at the true boundary of the resonator don't affect the symmetry pairing. For each symmetry 
choice one can evaluate the phase accumulated in the reduced resonator at each bounce, assigning a phase shift tt 
to each bounce off a "Dirichlet" internal boundary, and zero phase shift for each bounce off a "Neumann" internal 
boundary. If two symmetry types lead to the same final phase shift (modulo 27r) then those two symmetry types will 
be paired and quasi-degenerate, otherwise not. A subtle issue is the question of how to count bounces at the corner 
between two boundaries. The answer is that the semiclassical method really sums over orbits nearby the PO which 
will then hit both boundaries and experience the sum of the two phase shifts. 

We will illustrate this rule for the case of the bow-tie in the quadrupole. The symmetry reduced PO is shown 
in the last column of table O. It has one corner bounce, one bounce on the X axis and two boundary bounces. 
The boundary bounces don't matter as they will give the same phase shift for all symmetry types. The X axis 
bounce will give phase shift for the + symmetry of oz and ix for the — symmetry. The corner bounce sums the 
two shifts and gives: (+, +) — * 0, (+, — ) — * tt, (— , +) — * tt, (— , — ) — > 2-k. Adding these two shifts modulo 27r gives 
[(+, +), (+, — )] — > 0, [(— , — ), (— , +)] — > tt corresponding to the symmetry pairing we found above. In table ||, these 
two rules are applied to a number of relevant orbits in the quadrupole. It should be emphasized however that the 
group-theoretic projection method combined with the quantization rule which we illustrated in this section |IVB| will 
work for any symmetry group and the rules that we have stated are just useful shortcuts. 
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C. Evaluation of Mode Splittings 



The symmetry analysis above can only determine the existence of quasi-degenerate multiplets with small splittings, 
it cannot estimate the size of these splittings. In the phase space picture the splittings we are discussing come from 
tunneling between distinct periodic orbits, referred to as " dynamical tunneling" in the quantum chaos literature |55j . 
Techniques have been developed in that context for evaluating the splittings and we now apply those to the system 
we are considering. 

Dynamical tunneling in integrable systems is essentially similar to the textbook example of tunneling in one di- 
mension (note that conservative dynamics in ID is always integrable), and the splittings in such a case come from 
first order perturbation theory in the tunneling matrix element through an effective barrier, just as they do for the 
one-dimensional double well potential. Systems of the type we are considering, with mixed dynamics however show a 
very striking difference. It has been found by both numerical simulations and analytic arguments that the dynamical 
tunneling splittings in mixed systems are typically many order of magnitude larger than found for similar but inte- 
grable systems |36|, |37j (e.g in quadrupole vs. elliptical billiards). This difference can be traced to the mechanism of 
"chaos-assisted tunneling" (CAT) . As opposed to the "direct" processes when the particle (ray) "tunnels" directly 
from one orbit to the other, the CAT corresponds to the following three-step process: (i) tunneling from the periodic 
orbit to the nearest point of the chaotic "sea" , (ii) classical propagation in the chaotic portion of the phase space 
until the neighborhood of the other periodic orbit is reached, (iii) tunneling from the chaotic sea to the other periodic 
orbit. Note that the chaos-assisted processes are formally of higher order in the perturbation theory. However the 
corresponding matrix elements are much larger than those of the direct process. This can be understood intuitively 
as the tunneling from the periodic orbit to the chaotic sea typically involves a much smaller "violation" of classical 
mechanics and therefore has an exponentially larger amplitude. 




FIG. 7: Schematic indicating a direct tunneling process (black arrow) and a chaos-assisted tunneling process (yellow arrow) 
which would contribute to splitting of bow-tie doublets. 

The contribution of chaos-assisted tunneling can be evaluated both qualitatively and quantitatively using the so 
called "three-level model" |3(| where the chaotic energy levels (the eigenstates localized in the chaotic portion 
of the phase space) are represented by a single state Ec with known statistical properties. The straightforward 
diagonalization of the resulting matrix yields |3q| 



Ai?cAT 



\Vrc\ 2 
Er — Eq 



(58) 



where Er is the semiclassical energy of the "regular" states (localized at the periodic orbits) which does not include 
the tunneling contribution, and Vrc is the corresponding coupling matrix element with the chaotic state. The 
resonant denominator in Eq. ( |58| ) leads to strong fluctuations of the CAT-related doublets, as is found in numerical 
simulations |38| . The average behavior of the splittings however is determined by the matrix element Vrc- By virtue 
of the Wigner transformation J|6| it can be shown 39 that | Vrc 1 2 is proportional to the overlap of the Wigner 
transforms of the "regular" and "chaotic" states: 



\VrcV oc 



d sin x W c (</), sin \) W R ((/>, sin \) 



(59) 
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Assuming that, as required by Berry's conjecturepO[, on the average the Wigner function of a chaotic state is equally 
distributed across the chaotic portion of the phase space, and using the analytical expressions for the regular eigenstates 
calculated earlier, we find 



(A^cat) oc exp(-„4fci? ) 



(60) 



where A is the area in the Poincare Surface of Section (in (</>, sin x) coordinates) occupied by the stable island 
supporting the regular eigenstate. Note that Eq. (|o|) holds only on average, since chaos-assisted tunneling always 
leads to strong fluctuations of the splittings which are of the same order as the average |3fl] . 
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FIG. 8: The numerically determined splittings of bow-tie doublets for a closed quadrupole resonator with e = 0.14 (black dots) 
vs kR; the red line denotes the prediction of Eq. (^) for the average splitting, the blue line an estimate of the splitting based 
on the "direct" coupling. Note the large enhancement due to chaos-assisted tunneling and the large fluctuations around the 
mean splitting. 

In Fig. |§| we show the numerically calculated splittings for the bow-tie resonances in a resonant cavity with a 
fixed quadrupolar deformation e = 0.14, for different values of kR. Note that although there are large fluctuations 
in the numerical data (as previously noted), the data are consistent with Eq. (|60|), while a calculation based on the 
"direct" coupling severely underestimates the splittings. Unfortunately, due to the large fluctuations in the splittings, 
knowledge of the average splitting size does not accurately predict the splitting of a specific doublet. Note also that 
small violations of symmetry in the fabrication of the resonator may lead to much larger splittings then these tunnel 



splittings; such an effect was recently observed for triangle-based modes of GaAs ARC micro-lasers 19 



For the bow-tie orbit the quasi-degeneracy is associated with time-reversal symmetry, but as discussed above, it can 
equally well be associated with spatial symmetries and indeed this was the situation in the first work on dynamical 
tunneling |55j . Splittings of this type can also be adequately described using the framework developed in the present 
section. 
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V. CONCLUSIONS 

We have generalized gaussian optical resonator theory to describe the resonances associated with stable periodic 
orbits of arbitrary shaped two-dimensional dielectric resonators using the parabolic equation method. The corre- 
spondence to ray optics emerges naturally when imposing the boundary conditions at the dielectric interface which 
leads to the appropriate ABCD matrices for reflection, transmission and propagation. For a perfectly-reflecting 
cavity one gets quantized solutions at real values of k localized around the PO with mode spacings given by 
Afczong = 2ir/L, Aktrans = <p I L where L is the length of the PO and ip is the Floquet phase associated with the 
eigenvalues of the monodromy matrix (round-trip ABCD matrix). For a dielectric cavity one finds similarly localized 
quasi-bound solutions at quantized complex values of nk; in this case the imaginary part of nk is determined by 
the Fresnel refractive loss at each bounce of the PO. Within this approximation the mode spacings are unchanged 
from the closed case (except for the trivial factor of n). These regular modes coexist in a generic resonator with 
more complicated modes associated with the chaotic regions of phase space. Generalization of our results to the 
three-dimensional case appears straightforward for the scalar case and one expects only to have three-dimensional 
versions of the ABCD matrices enter the theory leading to some difference in details. More interesting would be the 
inclusion of the polarization degree of freedom, which seems possible in principle, but which we haven't explored as 
yet. 

We noted that for a cavity with discrete symmetries one will typically be able to construct several symmetry-related, 
nominally degenerate solutions of the wave equation for each PO. However group-theoretic arguments indicate that 
the these solutions cannot be exactly degenerate and lead us to construct symmetrized solutions which form quasi- 
degenerate multiplets. We presented a construction and then two simples rules for calculating the quasi-degeneracy of 
these multiplets and their symmetry quantum numbers. In the final section of the paper we show that the splittings 
of these multiplets are much larger than expected, due to the phenomenon of "chaos-assisted" tunneling, and estimate 
the average splitting in terms of classical quantities. 

There are several limitations of this work which we hope to address in future work. One obvious shortcoming is 
the prediction of zero width modes for the dielectric cavity if the underlying PO has all of its bounces above total 
internal reflection. Internal reflection is not perfect for these systems at any angle of incidence for two reasons. First, 
as these solutions describe gaussian beams with some momentum spread, every solution should have some plane- wave 
amplitude at an angle of incidence for which it can be partially transmitted. This type of correction exists even for a 
gaussian beam incident on an infinite planar surface and leads to an outgoing beam direction which can be significantly 
different from that predicted from SnelPs law; we have analyzed the beam deflection for this case recently [[ll]. For 
the case of the bow-tie modes of the dielectric resonator the effect of the momentum spread in inducing a finite width 
was also evaluated using a semiclassical method in ref. jl2). Second, due to the curvature of the interface in such 
resonators, there will always be some evanescent leakage, which we can think of as due to direct tunneling through the 
angular momentum barrier as is known to occur even for perfectly circular resonators. We have evaluated this type 
of correction also recently It is still not clear to what extent these effects can be accounted for systematically 
within a generalization of our approach here, e.g. to higher order in kl. One possibility we are exploring is that a 
generalized ray optics with non-specular effects included can describe the open resonator and its emission pattern. 
Both experiments ||, [l7j and numerical studies [[l7|, |4lJ demonstrate that for kl not too large (~ 50— 100) these higher- 
order effects must be taken into account. 

This work was partially supported by NSF grant DMR-0084501. 
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